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In magnetic fields applied parallel to the anisotropy axis, the relaxation of the magnetization 
of a Mni2-acetate single crystal measured for difi'erent sweep rates is shown to collapse onto a 
single scaled curve. The form of the scaling implies that the dominant symmetry-breaking process 
that gives rise to tunneling is a locally varying second-order anisotropy, forbidden by tetragonal 
symmetry in the perfect crystal, which gives rise to a broad distribution of tunnel splittings in a real 
crystal of Mni2-acetate. Different forms applied to even and odd-numbered steps provide a clear 
distinction between even step resonances (associated with crystal anisotropy) and odd resonances 
(which require a transverse component of magnetic field). 

PACS numbers:75.45.+j,75.50.Xx 



Single-molecule magnets are organic materials which 
contain a large (Avogadro's) number of identical 
magnetic molecules; ([Mni2 0i2(CH3COO)i6(H2 0)4]- 
2CH3COOH-4H20), generally referred to as Mni2- 
acetate, is a particularly interesting and much-studied 
example of this class. The Mni2 clusters are composed of 
twelve Mn atoms tightly coupled to give a sizable S* = 10 
spin magnetic moment that is stable at temperatures of 
the order of 10 K and below [Q. These identical weakly- 
interacting magnetic molecules are regularly arranged on 
a tetragonal crystal. 
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Magnetization versus longitudinal magnetic field 
for a Mni2 sample starting from a demagnetized state; data 
are shown for the same sweep rate at difi'erent temperatures 
ranging from 0.500K to 0.240K. The curves overlap indicat- 
ing that relaxation is from the ground state, as shown for the 
double- well potential illustrated in the inset. 



As illustrated by the double well potential shown in 
the inset to Fig. 0, strong uniaxial anisotropy (of the 
order of 65 K) yields a set of energy levels correspond- 
ing to different projections m = ±10, ±9, ,0 of the 

total spin along the easy c-axis of the crystal. Measure- 
ments below the blocking temperature of 3 K have 
revealed a series of steep steps in the curves of M ver- 
sus H at roughly equal intervals of magnetic field due to 
enhanced relaxation of the magnetization whenever lev- 
els on opposite sides of the anisotropy barrier coincide 
in energy. Below « 0.56K the magnetization curves are 
independent of temperature, and the tunneling proceeds 
from the ground state of the metastable well (see inset 
to Fig. 0). 

The spin Hamiltonian for Mni2 is given by: 



H 



-DSl~g,^XBH,S,-ASl 



(1) 



where D ~ 0.65 K is the longitudinal anisotropy, the sec- 
ond term is the Zeeman energy with gz ~ 1.94, and the 
third on the right-hand side represents the next higher- 
order term in longitudinal anisotropy. In order for tun- 
neling to occur, the Hamiltonian must also include terms 
that do not commute with Sz- In a perfect crystal, the 
lowest transverse anisotropy term allowed by the tetrag- 
onal symmetry of Mni2 is proportional to {S\ -\- St.)- 
For ground state tunneling, such a term only permits ev- 
ery fourth step. In contrast, all steps are observed with 
no clear differences in amplitude between them. This 
suggests that transverse internal magnetic fields, which 
would allow all steps to occur on an equal footing, pro- 
vide the dominant symmetry-breaking term that drives 
the tunneling in Mni2. However, dipolar fields |0-0| and 
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hyperfine interactions y,g|j9| are too weak to cause the 
rapid tunneling rates observed; the nature of the effects 
responsible for tunneling in Mni2 has remained an open 
question. 

In this paper we report data obtained for the relaxation 
of the magnetization of Mni2 in a swept field for differ- 
ent sweep rates. We show that a scaling form recently 
proposed by Garanin and Chudnovsky pO| , pl| , who con- 
sidered the effect of crystal dislocations, yields an ap- 
proximate collapse of all the data onto a single curve. 
The form of the scaling function corresponds to tunnel- 
ing due to second-order transverse anisotropy that varies 
throughout the Mni2 crystal with a very broad distribu- 
tion. Departure from perfect scaling is observed that is 
associated with a small admixture of tunneling due to 
other symmetry-breaking terms, presumably transverse 
internal magnetic fields. Thus 

n = + E{Sl - SD - H^S^. (2) 

with E ~ E{x,y,z) and Hx — Hx{x,y, z) varying from 
point to point in the Mni2 crystal. Our results imply that 
the dominant term responsible for tunneling is second 
order anisotropy which, although prohibited in a perfect 
crystal, is present and significant in real crystals of Mni2. 

The magnetization of small single crystals of Mni2- 
acetate was determined from measurements of the local 
magnetic induction at the sample surface using 10 x 10 
/im^ Hall sensors composed of a two-dimensional electron 
gas (2DEG) in a GaAs/AlGaAs heterostructure. The 
2DEG was aligned parallel to the external magnetic field, 
and the Hall sensor was used to detect the perpendicular 
component (only) of the magnetic field arising from the 
sample magnetization Jl2[ . 
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FIG. 2. Normalized magnetization curves at 0.24 K for 
sweep rates ranging from 1.30 x 10"^ T/s to 3.22 x 10"^ T/s. 
The plateaux indicate the cumulative fraction of molecules 
that have tunneled after an energy level crossing. For clarity 
only a partial data set is presented here. The inset illustrates 
the definitions of Tjv and -Rjv discussed in the text. 
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FIG. 3. The fraction of molecules, Rn, that remain in the 
metastable well following the A*'*'' level crossing plotted as a 
function of the logarithm of {dH^/dt)~^ . 

The magnetization of a Mni2-acetate crystal normal- 
ized by the saturation value is shown in Fig. ^ for dif- 
ferent sweep rates from 1.30 x IQ-^ T/s to 5.28 x 10"'' 
T/s. Each curve was obtained starting from a demagne- 
tized state by cooling the sample in zero magnetic field 
from above the blocking temperature to a low tempera- 
ture where pure ground state tunneling is observed. 

Since the magnetization has been normalized by its 
value at saturation, the data shown in Fig. ^ represents 
the cumulative fraction of molecules that have tunneled 
from the metastable into the stable well. The normal- 
ized magnetization at a plateau, labeled T^r in the in- 
set of Fig. 1^, thus represents the cumulative fraction of 
molecules that have tunneled (relaxed) following an en- 
ergy level crossing. By the same token, = 1 — T/v is 
the fraction of molecules that remain in the metastable 
well after the A^*'' level crossing. Fig. ^ shows Rn versus 
the logarithm of the inverse field sweep rate, {dHz/dt)~^ . 

For a single molecule, the probability of remaining in 
the metastable well is given by the Landau-Zener formula 
Pn = exp(-7rA^/2vjv) [|,|lO|,[l]jl||| , where An is the 
level splitting of the A^*'* resonance and is the energy 
sweep rate defined by = {gzfiQ'h/k'^){2S ~N)dHz/dt. 
To date, all attempts to account for the measured prob- 
abilities using the Landau-Zener method have required 
unreasonably large tunnel splittings. As we show below, 
a distribution of tunnel splittings, Atv.i (where i denotes 
the ith molecule), provides a more accurate description 
of tunneling in real crystals of Mni2-acetate. 

For a distribution of tunnel splittings, A^r^, the 
probability that a spin remain in the metastable well 
must be averaged over all the molecules: {PN,i) = 
^ exp(— ttA^ j/2u^), where Nt is the total number 

i 

of molecules. If the distribution is very broad, then (Pna) 
is best examined on a log scale, where an exponential 
looks like a step function, so that exp(— ttA^ j/2w^) « 
0(1 — ttA^ ■/2i'^). This means that for a fixed field 
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sweep rate, (dHz/dt), those molecules that have tunnel 
splittings obeying: 



(3) 



will remain in the metastable well for each N. In essence, 
Rn = {Pn,:) = ]^ E - 7rA%^j2v^). (4) 

i 

and each curve in Fig. || denotes the fraction of 
molecules that remain in the metastable well after each 
energy resonance because they tunnel too slowly p8| . 
Since Ajv.i is greater than Ajv~i,i, each step probes a 
set of molecules in the crystal belonging to a different 
range of the (initial) distribution of tunnel splittings. 

Eq. H shows that varying the sweep rate changes the 
fraction of molecules that remain in the metastable well 
after the field sweeps through a particular energy res- 
onance. In fact, Rn the fraction of molecules which 
satisfy Eq. up to the threshold condition ttA^ ^ — 2vj^. 
This suggests that varying the sweep rate provides a 
method for probing the distribution of tunnel splittings. 

Chudnovsky and Garanin |10 1^ have recently consid- 
ered a distribution of tunnel splitting due to transverse 
anisotropics caused by crystal dislocations. The formal- 
ism developed in references [^0|,^ is applicable regard- 
less of the physical origin of the distribution, provided it 
is logarithmically wide. Neglecting the quartic longitudi- 
nal anisotropy term AS^, the tunnel splitting of molecule 
i due primarily to second-order transverse anisotropy is 
given by: 



where 
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VN9N 



2D 



2D 



(25-Af)!(2S)! 
(TV)! 



(5) 



Transverse 



[(2S-Af-2)!!]2 

anisotropy generates selection rules that only allow even- 
numbered resonances. The leading contribution to even- 
numbered resonances involves ^ S — N/2 virtual 
transitions with Am — ±2; here 77^ — 1. However, odd 
numbered resonances could also occur due to the pres- 
ence of internal transverse magnetic fields of hyperfine 
or dipolar origin, or transverse fields due to the defects 
suggested by Chudnovsky and Garanin which locally tilt 
the easy axis so that the applied longitudinal field has a 
transverse component The largest contribution to 

odd- numbered resonances entails £,j^ = S — {N — l)/2 
virtual transitions with Am — ±2 due to the transverse 
anisotropy and a single virtual transition with Am = ±1 
due to the transverse field. One therefore expects that 
the tunnel splitting for even and odd N should be roughly 
comparable, as is observed. For odd N, the tunnel split- 
ting has the same form as Eq. ^. However, in this case, 
rjj^ = CN/2 and = 5 - (TV - l)/2, where C is an 
adjustable parameter of order 1. 

Rearranging Eq. ^ and employing the threshold con- 
dition, ttA^ j = 2i!jy, we find that 



In (I 1 /2D) = - In 
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N 



Hn^X (6) 



is independent of N . This implies that all the curves 
for i?7v shown in Fig. should collapse onto a sin- 
gle curve when plotted as a function of X (if second- 
order anisotropy is responsible for the tunneling). This 
is shown in Fig. IJ. 

We point out that scaling requires that different forms 
be used for even and odd- numbered steps, providing 
a clear distinction between even resonances (associated 
with crystal anisotropy) and odd resonances (which re- 
quire a transverse component of magnetic field) . For the 
odd resonance scaling a good fit was obtained for C = 1. 




FIG. 4. The fraction of molecules, 7?jv, that remain in the 
metastable well following the N^'^ level crossing plotted as 
a function of the scaling parameter X. The solid continu- 
ous curve is a best fit to the data using the error function. 
The inset illustrates the Gaussian distribution of transverse 
anisotropics determined by taking the derivative of this curve. 

Although the scaling obtained is of good quality, devi- 
ations should and do occur. This is due to the fact that 
the scaling function X was calculated exactly for even- 
numbered steps involving transverse anisotropy, while 
the expression for X is only approximate in the case 
of odd-numbered steps requiring an admixture of effects 
due to transverse magnetic fields. It should be noted 
that whenever more than one process contributes, per- 
fect scaling should not occur. 

The same formalism can be applied for the case of a 
distribution of tunnel splittings due instead to a trans- 
verse field that varies throughout the sample. In this case 
the tunnel splitting has the form: 



^'i,N = Sat 



'2D 



{2S-N) 



(7) 



If the tunnel 



whpyp n' - 2D / (25-jV)!(25)! 

wnere gj^ - [(^2S-N-iy.]'^ \J {Wy. • 
splitting were due to transverse field alone the data 
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should scale when plotted as a function of the scaling pa- 
rameter, - In {H,,^/2D) = - In {^[^mI^'n) /(25 - N) = 
X' . As can be seen in Fig. ||, the data do not scale by this 
procedure. We conclude that tunneling in Mni2-acetate 
is primarily due to transverse anisotropy, as evidenced by 
the approximate scaling shown in Fig. |j, with a small ad- 
mixture of tunneling due to transverse field which allows 
the odd numbered resonances and gives rise to deviations 
from perfect scaling in Fig. |[ 

Since T^v represents the cumulative fraction of 
molecules that have tunneled after the A^*'' and all pre- 
vious crossings, the negative derivative ol Rm with re- 
spect to X represents the distribution, /, of transverse 
anisotropy. The collapse of Rn onto one curve sup- 
ports the assumption that, rather than being the same 
throughout the crystal as is usually assumed, the tun- 
nel splittings vary locally within the Mni2 crystal with a 
very broad distribution. The fact that the universal re- 
laxation curve follows the error function as shown by the 
solid continuous curve in Fig. || indicates that the distri- 
bution is approximately Gaussian The distribution 
is shown in the inset of Fig. ^ 




FIG. 5. The scaling under the assumption that the tunnel 
splitting is due to transverse field alone. 

To summarize, the relaxation of the magnetization of 
Mni2 measured in longitudinal magnetic fields at differ- 
ent sweep rates collapses onto a single scaled curve. The 
form of the scaling implies that the dominant symmetry- 
breaking process that gives rise to tunneling is a locally 
varying second-order anisotropy, forbidden by tetrago- 
nal symmetry in the perfect Mni2-acetate crystal, which 
gives rise to a broad distribution of tunneling split- 
tings. Different forms applied to even and odd- numbered 
steps provide the first clear observation of a distinction 
between even step resonances (associated with crystal 



anisotropy) and odd resonances (which require a trans- 
verse component of magnetic field). 
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